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Abstract 

We establish a central limit theorem for partial sums of station- 
ary linear random fields with dependent innovations, and an invari- 
ance principle for anisotropic fractional Brownian sheets. Our result 
is a generalization of the invariance principle for fractional Brownian 
motions by Dedecker et al. 9 to high dimensions. A key ingredient 
of their argument, the martingale approximation, is replaced by an 
m-approxiniation argument. An important tool of our approach is a 
moment inequality for stationary random fields recently established by 
El Machkouri et al. [TC] . 
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1 Introduction 

Consider a stationary linear random field {Ci}jez<* with stationary mean-zero 
innovations {Xjjjg^d: 

= j;a,_iXi,j GZ^ (1) 

where {ai}i^^d are a collection of real numbers such that YlieZ'^ < oo. 
In particular, we are interested in the partial sum Sn '■= 

bnj = Eie{i,...,n}d ai-j and bn = (Ejezd ^nj)^^^- Then 5„ = Ejezd bnjXj. 
We establish sufficient conditions for the following two problems: 



(i) When do we have a central limit theorem 

Sn Ylij&Z'l^n^jXj 



K bn 



Ar(0, a' 
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(ii) When do we have an invariance principle 



and what is the hmiting process G? 

These two problems have a long history. In the one-dimensional case 
{d = 1), they have been extensively investigated and many results are 
known. For example, when {Xi}i^z are independent and identically dis- 
tributed (i.i.d.), Davydov in the seminal work [6] established an invariance 
principle for fractional Brownian motion (we will recall the definition be- 
low). See also Konstantopoulos and Sakhanenko [21] for a sharper condition. 
When {Xi}i(^z are stationary, the central limit theorem have been considered 
by for example Peligrad and Utev [351 ESj , Wu and Woodroofe [17j , Mer- 
levede and Peligrad [30], among others. The invariance principle has been 
investigated for fractionally integrated processes, an important class of our 
model ([1]) (see e.g. Wu and Shao [l6] and references therein for more results 
from the point of view of fractionally integrated processes), and recently by 
Dedecker et al. [9] in the general setting. 

However, very few of the corresponding results in high-dimensional case 
(d > 2) are known, with the notable exceptions of Surgailis |38) and La- 
vancier |28j. Surgailis [38] established general convergence results for (func- 
tionals of) linear random fields with independent innovations to self-similar 
random fields. Lavancier |28] investigated linear random fields with de- 
pendent innovations via a spectral convergence theorem (see also Lang et 
Soulier [25]). In [28], the limiting objects are described in terms of lin- 
ear mappings, and fractional Brownian sheets appear explicitly in the limit 
only when {Xijjg^d are i.i.d. random variables with zero mean and finite 
variance. 

Recah that, for a given vector H = {Hi, . . . ,Hd) G (0,1)*^, a (real- 
valued) fractional Brownian sheet = {B^jjgp,!]'' with Hurst index H 
is a real- valued mean-zero Gaussian random fields with covariance function 
given by 

E(Bf Bf ) = n ^ {sf^ + if ^ - K - t,\'''^), s, t G [0, 1]^ . (2) 

q=l 

When d = 1, is the fractional Brownian motion with Hurst index Hi. 
Fractional Brownian sheets play an important role in modeling anisotropic 
random fields with long-range dependence (see e.g. Doukhan et al. [14j and 
Lavancier [27]). They also arise in the study of stochastic partial differ- 
ential equations (e.g. Hu et al. [23J and 0ksendal and Zhang [E3]). The 
investigation of their sample path properties is another active research area 
(e.g. Xiao t48j). 
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When Hq = 1/2 for all q = in (l2|), the fractional Brownian 

sheet is just the (multiparameter) Brownian sheet, whence the former is a 
natural extension of the latter. Central limit theorems for stationary random 
fields and weak convergence to Brownian sheets have been considered by 
several authors, including Basu and Dorea [Tj, Bolthausen [3], Nahapetian 
and Petrosian [33] , Nahapetian [32] , Poghosyan and Roelly [37] , Dedecker [71 
|H], El Machkouri |15], Cheng and Ho ^5], Wang and Woodroofe [H] and 
El Machkouri et al. [16] . among others. In particular, people have recently 
found that in order to establish a central limit theorem for stationary random 
fields, a convenient way is to approximate the dependent random variables 
by m-dependent ones: the central limit theorems for m-dependent random 
variables have been known since Hoeffding and Robbins [22j, and it remains 
to show that the difference becomes arbitrarily small as m increases. This m- 
approximation method has been successfully applied to stationary random 
fields by Wang and Woodroofe [3T] and El Machkouri et al. [16], under 
different conditions measuring the dependence of random fields. At the 
same time, the m-approximation method has also been successful in one 
dimension. For example, recently Liu and Lin |29] applied this method 
to establish a strong approximation of stationary sequences by Brownian 
motions with optimal rates. 

This work develops an invariance principle for linear random fields with 
dependent innovations. Our result, as an invariance principle for fractional 
Brownian sheets, improves the ones by Surgailis [SS] and Lavancier [28] . as 
they both required innovations to be i.i.d. At the same time, their results 
are more general in the sense that they include weak convergence results 
with other limiting objects, while the ours does not. 

Our result can be seen as an extension of Dedecker et al. [2] to high 
dimensions. The main difference is that we replace their martingale approx- 
imation method, which seems difficult to be generalized to high dimensions, 
by the m-approximation method. Besides, to characterize the weak depen- 
dence of innovations, we apply the physical dependence measure introduced 
by Wu [43] and extended to random fields by El Machkouri et al. ^6]- In 
particular. El Machkouri et al. proved a moment inequality of weighted 
partial sums of {^ijjg^d, which is of significant importance in the analysis 
of random fields by m-approximation. 

Another crucial assumption for our results is a product structure for the 
coefficients: 

d 

fli = flj, {q) ,i eZ'^ , 

q=l 

and are square-summable real numbers for each q = 1, . . . ,d. In 

particular, the product structure allows us to extend the idea of coefficient- 
averaging by Peligrad and Utev [36] to high dimensions. It also plays an 
important role in the analysis of asymptotic covariance structure. 
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The product structure of the coefficients is an reasonable assumption and 
it was assumed in [28] to have an invariance principle for fractional Brownian 
sheets. To see that it is not a restrictive assumption, recall the product form 
of the covariance formula ([2]), and the fact that fractional Brownian sheets 
have the following stochastic integral representation with a product kernel 

where W = {Ws}^g]gd is a standard Brownian sheet, gHq{tq,Sg) = {{tg — 
Sq)^)^'>~^^'^ — ((— with r-|_ = max(r, 0) and hh is a normalizing 
constant (see e.g. [38]). 

At last, we point out that a more general problem is to consider partial 
sums of {i^(^i)}ieZ'* for some function ii' : M — M (we address the case 
K(x) = x), including the important case of empirical processes. A closely 
related problem is to consider the asymptotic properties of partial sums of 
{ii'(Zj)}jg2d where {Zi}j^^id is a Gaussian random field. Such results have 
significant impacts in statistics theory. These problems again have been ex- 
tensively investigated in one-dimensional case (see e.g. Taqqu [39], Giraitis 
and Surgailis [18,19], Ho and Hsing [201 [2l] , Wu [121 gl] , to mention a few), 
while high-dimensional cases have been much less considered (see e.g. Do- 
brushin [11] and Dobrushin and Major [12] . Surgailis [38], Doukhan et al. |13) 
and Lavancier [27]). To extend our results to general functionals of linear 
random fields is beyond the scope of this paper, but worth investigating in 
the future. 

The paper is organized as follows. We will provide the background on 
the physical dependence measure in Section [2l A central limit theorem 
(Theorem and an invariance principle (Theorem [3]) are established in 
Sections [31 and H) respectively. Discussions on related works are provided in 
Section [5l 

2 Preliminaries on physical dependence measure 

Consider stationary random fields {Xi}^^^d of the following form 

Xi=g{ei.j -.j (3) 

where g : M^'' — )• M is a measurable function and {ej}jgz<* i.i.d. random 
variables. Throughout this paper, we assume that KXq = 0. 

El Machkouri et al. [I6j suggested to measure the dependence of {Xi}^^j^d 
as follows. Let e* = {e*}i^id be a random field coupled with e, defined by 
e* = ej for all i 1/ \ {0} and being a copy of eg independent of e. 



4 



Set X* = g{e*_j : j € Z*^) and define the physical dependence measure of 

A,^Ap{X)=Y,\\X^-X*\\^. (4) 

El Machkouri et al. derived a central limit theorem and an invariance prin- 
ciple under the condition that Ap < oo for certain p > 2. In particular, the 
following result is useful for our purpose. 

Theorem 1 (El Machkouri et al. |16]). (i) Let {ai}^^^d be a family of real 
numbers. Then for any p > 2, 

\\^aiXi\\ < (2p J^a?) Ap. (5) 

(a) A2 < 00 implies that ^f^^^d |EXoXfc| < 00. 

In fact, dS]) was proved for {ai}j^^^d having finite non-zero numbers, but the 
extension is immediate. 

Our results require stationary random fields {Xi}^^^d to satisfy Ap < 00 
for some p >2. El Machkouri et al. flE] provided several such examples. 

Example 1. Consider {^j}jgz<* form of functional of linear random fields: 

where 5 is a Lipschitz continuous function and the coefficients {V'ijigZ'' sat- 
isfy X^jggd IV'jl < CO. If eo € LP for p > 2, then Ap < 00 (p^, Example 1). 
Note that this class of random fields include linear random fields. Another 
class of non-linear random fields are the Volterra fields ( [16] , Example 2) . 

Remark 1. Moment inequalities play an important role in establishing 
asymptotic results for random fields. Dedecker [8] also established a similar 
moment inequality for random fields, under a different condition of weak 
dependence. In principle, in order to establish asymptotic normality one 
should expect to control, for finite subset T C Z'^, || Ylier -^iWp — 
for some p > 2. Wang and Woodroofe [41] established such an inequality 
for r in form of rectangles, under a different condition of weak dependence. 

Remark 2. In the literature of (one-dimensional) stationary sequences, such 
a condition Ap < 00 on the weak dependence is often referred to as of pro- 
jective type. An advantage of projective- type conditions is that, they often 
lead to easy-to- verify conditions for the asymptotic normality of various sta- 
tionary processes arising from statistics and econometrics. For more on the 
projective- type conditions in one dimension, see for example Wu [431145] and 
Merlevede et al. [31j. For other types of conditions on weak dependence, see 
for example Bradley [1] and Dedecker et al. [10] . 
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3 A central limit theorem 



We first establish a central limit theorem for triangular array in form of 

Sn = Sn{X,b):=Y,bn,jXj, (6) 

with general regw/ar coefficients {bn,j}n,j to be defined below. For each / S N, 
set rectangle blocks of size l'^ by 

Ik = Ik{l) = {ieZ'^:ige {Ikg + l,...,lkg + l},q = l,...,d} CZ'^, 

and define ^ 

^rr,k = JdYl ^" J ,n£'N,k£Z'^ (7) 

and Cn = {YlkeZ'i '^n k)^^'^- introduce the following definition in the spirit 
of the coefficient-averaging idea of Peligrad and Utev 



Definition 1. We say the coefficients {{bn,j}j(zzd : n G N} are regular, if 
6^ — > cxD and for each Z € N, 

i™o5^E E(V.-Cn,.)' = (8) 

hm ^ Yl Kj-'^lkl = (9) 
lim sup = 0. (10) 

Theorem 2. Consider Sn as in ([6]) with some regular coefficients {hn,j}n,j- 
//A2 < 00, then, cj^ := ^^^^dlEXoXfc < 00 and 

^^M{0,a'). (11) 

Theorem [2] leads to an answer to our first question. In particular, return 
to our problem of linear random fields with bnj = X^ieji n}<^ cii-j- It suf- 
fices to show such {bn,j}n,j are regular. Recall that the coefficients {aijig^d 
are assumed to have the product structure: 

d 

ai = Y[ai^{q), for some {ai^(g)}i^6z,g = 1, ••• ,d, (12) 

q=l 

with J^jez-i '^i < This suffices to establish the regularity of {bn,j}n,j- 

Corollary 1. If bnj = J2ie{i nY^i-j '"^^^^ {^i}ieZ'' having the product 
form ()12p . then {&n,j}n,j is regular. As a consequence, if in addition A2 < 
OO; then Sn{X,b)/bn =^ A/'(0,cj^) with cj^ = Ylj&'^'^i^Q^j) < 
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We first prove Theorem[2]through a series of approximations of S'n in (l6|). 
Tlie main tool is Theorem [TJ 

(i): m- approximation. First, fix m and let X = {Xi}^^^d (depending on m) 
denote the stationary random field obtained by 

Xi = E{Xi I ) , 

where J^™ = cr{€i-j : j G {— \m/2\ , . . . , [m/2j}'^). The so-obtained random 
field {Xjjjg^d is (m + l)-dependent, that is, for all i,j € Z*^, Xi and Xj are 
independent, if maXq=i^ \iq — jq\ > m + 1. 

Lemma 1. // Ap < oo for some p >2, then 

r \\SniX,b)-Sn(X,b% 

lim sup ; = . 

m^oo n On 

Proof. By Proposition 3 in |16] . 



\\Sn{X,b) - Sn{X,b)\\<(2p 6^,,y^'Af ) , 

where Ap™"^ is the physical dependence measure for X — X: 

AM = V -Z,) - (X, -Z,ri|p. (13) 
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By Lemma 2 in [16j, Ap < oo implies that limm— !>oo 

Af ) =0. □ 

In the sequel, let Ap = ^^g^d —J(*\\p denote the physical de- 
pendence measure of X. Observe that Ap < oo implies that Ap = 
Eie{-m,m}rf ll-^j - ^*j\\p < oo> for all m € N. 

(a) Coefficient-averaging. This procedure was introduced by Peligrad and 
Utev [36] in the one-dimensional case. For each / G N, recall the definition 
of Cn,k in dZ]). Set 6„j = Kjil) = Y^kez^ Cn,fcl{je/fc}- 

Lemma 2. For each m G N, ^ € N, if Ap < oo for some p>2, then 

^.^ \\Sn(X,b)-Sn(XM\v ^ Q 

n— s>oo 5n 

Proof. Apply Theorem □ and ([8]) to X. □ 

(^Mzj Big/small blockings. Define, for k G Z'^, Z € N, / > m + 1, 

4 = 4(0 = {i G Z^ : € {Ikq + 1, . . . , //c^ + / - (m+ 1)}, g = 1, . . . , d} C 4 . 

Set Yfc = Sjg7^ , G Z'^. Since {Xjjjg^d are (m + l)-dependent, by 
the construction of {Ik\k&i<i-, {^A:}fcez<* are i.i.d. random variables. Consider 
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Lemma 3. For each m € N, i/ Ap < oo for some p >2 and ([9]) holds, then 

Z— >oo n—^oo bji 

Proof. Observe that Sn{X ,b) - SniY, c) = J^kez^ Cn,k I]jg4\7j^ ^j- % The- 
orem (TJ 



\\Sn{X,b)-Sn{Y,c)\\p < {2p[/'^-(Z-(m + l))'^] J^c^ ,} 

= {(2p)[l-(i^(!^)^](,V„)}"^A,. 

Note that Q imphes hm„_j.oo I'^c'^/bn = 1, whence the desired result follows. 

□ 

(iv) Triangular array of weighted i.i.d. random variables. Now we establish 
a central limit theorem for SniY, c). Recall that Y depends on /, m G N. 



Lemma 4. For each mGN,l>m+l, if Ap < oo for some p >2 and (|T0 
holds, then 

with 

<i= E n (i - 

ie{m+l~l,...,l~m-l}<' r=l 

Proof. It is equivalent to prove a central limit theorem for X^jg^d CnjYj/cn 

with Yj = Yj/l'^/'^. By straight-forward calculation, 'E(y\)/1'^ = cr'^i. 
Then, (fTOjl yields the desired result. □ 



Proof of Theorem [H Combining Lemmas HHH yields that 

\\Sn{X,b)-Sn{Y,c)h „ 

lim iim iim sup = L) . 

m^oo l^co n^oo On 

By [in], proof of Theorem 3.1 therein, = limm^oo limz-5>oo /• The 
desired result follows. □ 

At last, we prove Corollary [TJ By the product form ()12p . we can write 

d d 
bn,j = Whn,j^{q) and 6n = (14) 

q=l q=l 
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with bnj,{q) ■■= Er=i«i-i,(9) and := [T.j^eAjS^)V''^ ' Accordingly, 

write c„,fc^(g) := Ylfj=k}\iKiM) ^nd c„(g) := 



for each l,k G 



[E 



kq& ^n,kq 



In the sequel, for the sake of simplicity, for j,k £ write bnj{q) = 
bn,jq{q) and c„^fc(g) = Cn^kqio) and view them as functions of jq,kq G Z, 
respectively. 

Proof of CorollaryUl Fix Z G N and recall that Cn^k depends on /. We first 
show supjg^d \cn,j\/cn — >■ as n — >■ GO. Write 



sup — 



< 



n 



sup 



\Cn,j{q)\ 
Cn{q) 



(15) 



and it converges to zero by the fact that sup^g^d \cnj{q)\/cn ^ as n ^ oo. 
It suffices to show sup^g^d \bnj{q)\/bn{q) — ?> as n — )■ oo for all g = 1, . . . , d, 
which was proved in Peligrad and Utev [35J, p. 448. 

Next, we show that the product form ()12p implies that dS]) and ([9]) hold. 
This result is an extension of Peligrad and Utev [36], Lemma A.l. We prove 
by induction. Note that now ([9]) becomes 



^ d d 



9=1 



0. 



(16) 



When d = 1, it was shown in [36j that ([5]) and (fTUj) holds for Ylkez < oo- 
Suppose dH]) and (fT6|) have been proved for d — 1. We prove ([8]) for d and 
the proof of (|16|) is similar and omitted. By the inequality that for any real 
numbers Og, f3q, q = 1, . . . , d. 



d-1 



d-1 



d-1 



5=1 q=l g=l 5=1 g=l 

we bound Efc EJnLi fen,,(9) - nil Cn,fc(g)]' < 2($i'^ + $i'^) with 



E E [ n - n cn,fc (^z)] ' X 62 (d) 



ki,...,ka_ljl,...,jd-i 9=1 



9=1 



and 



^i'^ = E E n<fc('?)EEK(^)-^"-'^'(^) 

A;i,...,A:fj_i ji,...,jd-i 9=1 ^d id 

d-1 

= ^'-'n^n(9)EEK^(^)-^".'^('^) 
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By induction, 



(1) 



E 



Y{t=\cn,k{q) 



62 



and similarly, ^>i^V&n ~ Efc, EjJ^n,i(rf) 
with d = 1. We have thus obtained 



Cn,kmVbl{d) 
for all d G N. 



= o(l), 
o(l) by 



□ 



4 An invariance principle 

We consider weak convergence in the space -D[0, l]'^ consisting of functions 
'continuous from above with limits from below' (see Bickel and Wichura [2] 
for details). For t G [0, 1]'^, consider Sn{t) = • • • Y}C=} ^^ ^ D[0, 1]'^. 

This time we have 

[nti] [nta] 

Sn{t) = ^ bnt,jXj with bnt,j ■= XI ' ' ' X] "'^'^ (^'^^ 
jeZd ii=l id='^ 

Theorem 3. Suppose there exists H G (0, 1)*^ such that 

lim ^^^^^ = s^^', for all s G [0,l],g = l,...,d, (18) 
and there exists p such that 

p > 2, p > max — and Ap < oo . (19) 

q=l,...,d Hq 

Then, {<S'n(t)/6ri}tg[o,i]d converges weakly in D[0, to the fractional Brow- 
nian sheet with Hurst index H . 

Remark 3. When (fT8|) holds for some H G (0, 1)*^ with maxg=i^,,,^rf i?"^ < 2, 
condition (|19|) becomes A2 < 00. Otherwise, we need to assume finite 
higher-than-second-order moment to establish the tightness. A similar phe- 
nomena was observed in the one-dimensional case ([S], Theorem 3.2). 

Example 2. Due to the product structure, it suffices to provide examples 
of {aiq}j,6Z such that (fT8]) holds for each q = 1,. . . ,d. Several examples 
have been provided in Dedecker et al. 0|, Examples 1-4. We summarize 
them below. 

(i) Fix a G (0,1/2), and set oq = l,ai = r(z + a)/(r(a)r(i + 1)) fori > 1. 
Then H = a + 1/2. 

(ii) Fix a G (0,1/2), and set m = {i + 1)"° - i"" for i > 1. Then 
H = 1/2 -a. 
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(iii) Fix a S (1/2, 1) and set ai ^ i for i > 1 with any slowly varying 
function / at infinity. Then H = 3/2 — a. 

(iv) Fix a > 1/2 and set ~ r^/'^{logiy°' for i > 1. Then H = 1. 

Example 3 (Fractionally integrated random fields). Case (i) in Example [2] 
above corresponds to the fractionally integrated random fields, generated 
by back-shift operators. Let Bg denote the back-shift operator on the q-th 
coordinate of the random fields {Xjjjg^d: BgXi = 
Then, for G (0,1/2), 

oo 

(/ - Bg)-"^X, := a.(9)X,,,...,,,_,,,-,_,,,^^„...,,, , j G Z^ 

with ai{q) defined in Example [21 (i). Thus, fractionally integrated random 
fields defined by 

= {I- Bi)-"^ •••(/- BdY'^'^Xj , j G 

fit in our model ([1]) with coefficients Oj = Y^q=iaiq{q)-,i G Z"^, where 
{^iq{<i)}i<^i+ corresponds to Bq as above and ai^{q) = for iq < 0. This gen- 
eralizes the fractional autoregressive integrated moving average (FARIMA) 
processes (see e.g. [16] and references therein) to random fields ([27J). Note 
that Wu and Shao [l6] also established an invariance principle for the so- 
called Type II fractional integrated processes, which are slightly different 
from our model ([1]). 

Theorem [3] follows as usual from the convergence of finite-dimensional 
distributions and tightness (see e.g. Bickel and Wichura [2]), which are 
proved below separately. 

Proposition 1 (Convergence of finite-dimensional distributions). Suppose 
A2 < 00 and for some H G (0, l)'^ (jlSp holds, then the finite- dimensional 
distributions of {<S'n(i)}te[o,i]'* converge to that of a fractional Brownian 
sheets {^t^}te[o,i]'' 'ujith Hurst index H . 

Proof. We start with some new notations following ()17p . This time the 
product structure of the coefficients yields, for n = (n(l), . . . , n((i)) G Z^J., 

d n{q) 

Kj = n ^n{q),j, (q) with bn(q)J^ (q) = fflj-j, (?) , j G Z'^ . (20) 
q=l 1=1 

Define accordingly 6„(g) = [Ej^g^ ^^(g) j_^(g)]^/^- Similarly define Cn,k{q) = 

T.f'=k}+iK{q),j,{'i)l^ and Cn,k = nq=ic„{g),fc,(9). As before, for n G Zf 
and j,k G Z'^, write 

Kjiq) = K(q),jg{q) and Cn,k{q) = Cn(q),k,{(l) , (21) 
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for the sake of simplicity. 

Fix m G N. Take arbitrary t^^\ . . . ,t^'^^ G [0,1]'^ and write = 
( lnt^{^\ [nt'-J^l ) G for r = 1, . . . , m. Then we can write 

m m 
r=l j(zid r=l 

By Theorem [21 it suffices to show that {hn,j = YlT=i Kbnr,j}n,j are regular 
and 



bn--= E {Y.^rbn.'^) ~ &^Var [ J] A.B^ (t^ ) 



(22) 



r=l 



r=l 



We first prove ([22]) . Observe that 

mm m m 

jeZ<* ^■=1 s=l r=l s=l j^ld 



ns,3 ■ 



For fixed r,s, recall that for j G Z^, by our notation ([2T]). bn^j{q) depends 
only on jq G Z. Then, 



d 



q=l 



We have thus shown that 



d 



mm ^ -I 



r=l s=l g=l 
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On the other hand, by the covariance formula ([2]), 

m 
r=l 

7TL TTL d 

r=l s=l q=l 

Now ([22]) follows from ([IHD by recalling that = nq=i bl{q). 

Next we check that {bnj}j^^d are regular. Accordingly define c^^k — 

^jeh^nj/l"^ and c„ = EfceZ'* c^,fc]^^^- Observe that Cn,k = E^i^rCn.j- 
Then, conditions ([5]), and ^ become 

^ ^ ^ (bn,j - Cn,kf ^ and i ^ ^ |6^j - c^j,| ^ as n ^ cx) . 

(23) 

The first part follows from the observation that for each r, {b-nj — Cn,k)'^ < 
^117=1 >^ribnr,j-Cnr,kf and EfceZ'' T^jeiS^nrj-Cnr^k)"^ /bl ^ as n ^ cx). 
To show the second part of ()23p . observe that 

l^nj ~ - X/ ^ l^"'-?' ~ '^ri,fc||^nj + C„,fc| 

< ( 5^ 5^ IVi - Cn.fcP) ^ ( 5^ 5^ IVi + C„,fcp) ^ , 

where the first term in the last product is of order o(6„) (by the first part 
of ([23]) ). while the second term 0{bn), whence (f23l) follows. 

At last, condition ([TO]) becomes lim„_s.oo sup^g^d |c„j |/c„ = 0. To see 
this, observe that the second part of (|23|) implies that I'^c^ ~ ^n- then 
follows from ([22]) and ([9]) that ~ Cc^ for some constant C > 0. The rest 
of the proof is similar to the control of p5|) and omitted. We have proved 
the regularity of {bnj}n,j and thus the proposition. 

□ 

Proposition 2 (Tightness). If there exists p such that ()19p holds, then the 
process {Sn{t) /bn}t^[Q ^d is tight in D[0,l\'^ . 

Proof. We will apply Lavancier [26j, Corollary 3. By slightly modifying the 
argument therein, it suffices to show that there exists constants /3 > l,p > 
0, C > 0, such that for all t € (0, 1)°' and n large enough, 

d 

\\Snmi<CKU*^,- (24) 

q=l 
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By Theorem [H for p > 2, 

d 



\Sn{t% = II ^nt,iX^ ^ (2p)^/'6ntA, = {2p)'l^ J] ^'1 ■ (25) 



Observe that (fT8|) means that hn{q) is regularly varying. Now, by Taqqu 
Lemma 4.1, for all 7^ > 0, there exists Cg > such that hnt/bn < Cqt^''~'^'', 
uniformly on [0, 1] for n larger than some Ug. Now, ()25p can be controlled 
by, for n large enough and some constant C, 

9=1 

If p satisfies (|19|) . then one can choose 7g > small enough so that p{Hq — 
7g) > 1. It then follows that dMl) holds with (3 = m.inqp{Hq - 7^) > 1. We 
have thus proved the tightness. □ 



5 Discussions 



We compare our results with Surgailis [38] and Lavancier [2H]- Surgailis 
proved more general results in the sense that he considered general func- 
tional of linear random fields, with independent innovations. However, it 
assumes finite moments of any order of the innovations, which is not neces- 
sary in the special case of invariance principle for fractional Brownian sheets. 
It is an interesting problem that whether the arguments shown here can be 
extended to the more general settings of [38j . 

Next, we compare our results and Lavancier [28]. The latter assumed 
the following assumption on the stationary random fields: 

HI The random field {Xi}^^^d is weakly stationary (i.e., with shift-invariant 
covariance structure), centered and has a bounded spectral density fx, and 
there exists a random field {^t}te{o,oo)''- such that 

[nil J [ntd] 

jl=l jli=l 

By TheoremlU (ii), A2 < cxo implies that the spectral density of {Xi}i^^d 
is bounded (see e.g. Fan and Yao [17J, Theorem 2.11 for the case d = 1). 
By [16j, Proposition 4, A2 < 00 also implies (|26p . Therefore, our assump- 
tion on the random fields is stronger than HI. Results in |28j are also more 
general in the sense that they do not assume the product structure of coeffi- 
cients, and they cover various limiting objects including fractional Brownian 
sheets. 
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However, as an invariance principle for fractional Brownian sheets, our 
setting is more general than [28]. Indeed, the limiting objects in |28] are 
described through a linear mapping (from L^{R'^) to L^{n,B,F)), which 
is not explicit in general ([28j, Remark 4). In particular, the fractional 
Brownian sheets can be interpreted from this characterization, only when 
{Xjjjg^d are i.i.d. and the coefficients are of product form. 

We conclude by comparing our assumptions on the coefficients with the 
ones in |28) , Theorem 5 (where strong white noise and product form of coef- 
ficients were assumed). Due to the nature of the spectral analysis approach 
therein, for each q, the coefficients {ai{q)}i^z were assumed to be the Fourier 
coefficients of certain function a^'^^ € L-^([— vr, vr]): 

a^'^\u) = Y,a,{q)e-^'^'^,q = l,...,d. 

Focus on a^'^\ujq) and omit the index q from now on. In [2^, Theorem 5 
(see also Remark 6), it was assumed that 

a{u) ~ C|cj|~", for some a £ (0, 1/2), C> . (27) 

By results on trigonometric series (see e.g. Zygmund [49j, Chapter V, The- 
orems 2.6 and 2.24), (f27|l is equivalent to assume aj ~ Cij"~^ as j — )• oo for 
some constant Ci > 0, which is a special case in Example [21 (iii). The other 
cases of Example [2] are not covered by ()27p . Therefore, our assumptions on 
the coefficient are more general than ^28j in the case of invariance principles 
for factional Brownian sheets. 
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